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Abstract 

A search for day-night variations in the solar neutrino flux resulting from 
neutrino oscillations has been carried out using the 504 day sample of solar 
neutrino data obtained at Super-Kamiokande. The absence of a significant 
day-night variation has set an absolute flux independent exclusion region in 
the two neutrino oscillation parameter space. 

14.60.Pq,26.65.+t,96.40.Tv,95.85.Ry 
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As a real time solar neutrino experiment, Super-Kamiokande can perform a wide range 
of time modulation studies of the solar neutrino flux. One motivation for these types of 
studies is an investigation of neutrino oscillation hypotheses. All solar neutrino observations 
]l|-[5|] have reported significantly lower fluxes than the expectations of standard solar models 
(SSMs) This difference is commonly referred to as the solar neutrino problem. Given 

the support of recent helioseismological observations |J, these SSMs look well-established 
and reliable. The difference between observations and predictions suggests some neutrino 
properties beyond the standard model of elementary particles. The most popular solution 
to the solar neutrino problem is neutrino oscillations, aided by matter enhanced oscillations 
in the Sun |Kj. In some regions of the parameter space for neutrino oscillations, matter 
enhanced oscillations within the Earth can lead to a regeneration of the measured neutrino 
flux passing through the Earth. This regeneration would produce a higher flux measured 
during night time relative to day time measurements. If such a day-night variation was 
observed, it would be strong evidence for neutrino oscillations. The amplitude of the day- 
night flux variation would determine the neutrino oscillation parameters, independent of the 
absolute flux uncertainties of the SSMs. 

Super-Kamiokande (SK) started taking data in April, 1996. SK has already confirmed 
the deficit of solar neutrinos ||. In this report, the total live time is increased to 503.8 
days (May 31 st , 1996 through March 25 th , 1998), and the total number of solar neutrino 
events found coming from the Sun is now 6823ti3o events above a threshold of 6.5 MeV 
in total energy of the recoil electron. Day-night (DN) variations are investigated with this 
high-statistics solar neutrino data sample and an updated flux value is presented. 

The SK detector is located at the Kamioka Observatory, Institute for Cosmic Ray Re- 
search, the University of Tokyo, in Gifu Prefecture, Japan, 137.32 degrees East longitude 
and 36.43 degrees North latitude. Due to the latitude, the nadir of the Sun can range be- 
tween ±0.974 in cosine (12.98 to 167.02 degrees) at the SK site. In this analysis, nadir {6 Z ) 
is defined as the angle between the negative z-axis of the detector coordinate system and 
the direction to the Sun (solar neutrino direction), where the cosine of the nadir is positive 
when the Sun is below the horizon. Solar neutrinos will penetrate different regions of the 
Earth, that are related to the nadir of the Sun by the following: mantle (0 < cos 6^ < 0.838), 
outer core (0.838 < cos# z < 0.981) and inner core (0.981 < cos #2 < 1). SK never sees solar 
neutrinos that pass through inner core of the Earth. However, the fact that the density 
of the outer core (10 ~ 12g/cm 3 ) is about double that of the mantle (< 5.5g/cm 3 ) may 
enhance the neutrino regeneration efficiency of the Earth. To look for effects associated 
with the core, the night sample was divided into 5 data sets according to nadir of the Sun 
at the time of the neutrino event, Nl (0 < cos6» 2 < 0.2), N2 (0.2 < cos6> 2 < 0.4) ... , and N5 
(0.8 < cos # 2 < 1.0). In the cases of Nl to N4, neutrinos pass only through the mantle of 
the Earth. N5 data contains the outer-core-penetrating neutrino sample. In the N5 period, 
the Sun has spent 80% of its time in the outer core region. 

The absolute energy scale of the SK detector has been calibrated precisely using a linear 
accelerator system (LINAC) at the detector O]. The LINAC injects electrons downward 
in direction, at several fixed points in the volume, and at certain times. Uniformity of the 
detector response in each nadir direction is important, since the data is divided according 
to the nadir angle. The uniformity in azimuthal angle is also important since the nadir 
and azimuthal angles of the Sun are correlated. In addition, N5 data can only be taken 
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during the winter time, so long term stability of the detector must be monitored. The 
precise LINAC calibration has to be extrapolated to the entire volume, in all directions, 
and at all times. Decay of spallation products induced by cosmic ray muons are used for 
this purpose. Spallation events accumulate in the data as a major background in the solar 
neutrino measurements (~ 600 events/day above 6.5 MeV threshold in a 22.5 ktons fiducial 
volume). They are distributed uniformly in volume, direction, and time, just as solar neu- 
trinos. Their beta decay energy spectra are distributed in the relevant energy region for 8 B 
neutrinos, making this background a good calibration and stability monitor. The spallation 
event sample was divided into subsets in time and direction (nadir and azimuthal angle) 
and the relative energy difference was found by comparing the spectral shape of a subset 
with that of the whole spallation sample as a reference spectrum. The magnitude of spread 
of the relative energy differencies is consistent with the expected statistical distribution of 
subsets and no systematic biases are seen in the test variables. A conservative systematic 
error for relative energy scale for each Nl to N5 data set is 0.5%. Due to the steepness of 
the recoil electron spectrum near the analysis threshold of 6.5 MeV, this small scale error 
is amplified and the relative flux error of each data subset (N1-N5) in the uniformity of the 
energy scale is estimated to be Since the statistical errors of the flux values are much 

larger than this scale error and the spread of the relative energy scales was consistent with 
a statistical distribution, we treated this scale error as uncorrelated among subsets in the 
following analysis. 

In order to be independent of the absolute flux values of 8 B solar neutrinos in SSMs, only 
the relative difference of DN flux values were used in the oscillation analysis. In this case, 
many of the systematic errors involving detector response cancel. The remaining systematic 
errors, in addition to the relative energy scale error, are estimated to be: ±0.1% (data 
reduction), ±0.4% (background subtraction), ±0.1% (live time calculation), and < 0.1% 
(trigger efficiency) for each night subset. 

Figure [I] shows the measured DN fluxes from 503.8 days of SK data. The flux is nor- 
malized to a value at one astronomical unit (AU) by a correction for the eccentricity of the 
Earth's orbit. Numerical values of the corrected fluxes are listed in Table [J] with relative sys- 
tematic errors. The total flux updated for 503.8 days of data is also listed in the table with 
an absolute systematic error same as previously presented ||. Obtained DN asymmetries 
are as follows: 



No significant DN variation nor N5 excess is seen in the data. 

The impact of these results was investigated within a two neutrino oscillation hypothesis 
for v e — > or u e — > v T (active neutrinos), and v e — > v s (sterile neutrino). A flux independent 
analysis was performed by treating the flux normalization factor a as a free parameter in 
the x 2 definition. 



N 



1 = 0.047 ± 0.042(stat) ± 0.008(syst) 



D 

N5 



0.055 ± 0.063(stat) ± 0.013(syst) 



(D, Nl--- N4) 
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where 4>i is the measured flux, 0° sci (sin 2 29, Am 2 ) is the effective flux for a given set of oscil- 
lation parameters derived from the ratio of the expected number of events with and without 
oscillations, cr, and u sys ^ are the statistical and systematic errors of i t/l -bin listed in table |. 
The x 2 value for the case of no oscillations is 7.4 with 5 degrees of freedom (dof), which 
corresponds to a 19% probability. The expected solar neutrino flux nadir dependencies for 
a set of oscillation parameters were obtained by a numerical calculation using models for 
the neutrino production point M, electron density in the Sun H, electron density in the 
Earth (PREM [JP2]] ) and 8 B neutrino spectrum [[HJ. Production points were integrated on a 



0.01R sun grid in the plane which contains the Sun and the Earth. The electron density dis- 
tribution in the Earth was calculated with charge-to-mass ratios (Z/A) of 0.468 for the core 



and 0.497 for the mantle [fL4|l . Neutrino trajectories in the Earth were integrated over 1000 
directions with 0.001 cos9 z steps, weighted by SK live time. Assuming neutrino incoherence 
at the Earth, the electron neutrino survival probability at the detector, Psei was obtained 
from independent calculations of Pi j2 (Am 2 /P, sin 2 29, r ) and Pi ej2e (Am 2 /P, sin 2 29, cos 9 Z ) 
using 

P S E = Pi Pie + P 2 P 2e = (1 - P 2 )(l " P 2 e) + P 2 P 2 e 

where Pi,P 2 are the probabilities to be v\,V2 at the surface of the Sun and ro is the pro- 
duction point of neutrinos in the Sun. P\ e , P 2e are the probabilities to be detected as a u e 
at SK if the neutrinos arrive as v\, z/ 2 , taking into account any possible regeneration in the 
Earth. Representative flux variations obtained from these calculations (for active neutri- 
nos) are shown in figure |l| for a typical set of parameters at a large mixing angle solution 
(sin 2 29 = 0.56, Am 2 = 1.2 x 10~ 5 e\^ 2 ) and for a set at a typical small mixing angle solution 
(sin 2 29 = 0.01, Am 2 = 6.3 x 10- 6 eV 2 ). 

The expected solar neutrino DN flux was calculated for points in the parameter space 
(10 -4 < sin 2 29 < 1, 10~ 8 < Am 2 < 10~ 3 eV 2 ) and compared to the measured values. 
Minimum \ 2 values of 5.2 were found at (sin 2 29 = 3.2 x 10™ 2 , Am 2 = 1.6 x 10~ 6 eV 2 ) for the 
v e —* case and 4.8 at (sin 2 29 = 0.25, Am 2 = 1.8 x lO^eV^ 2 ) for the v e — > v s case. These 
sets of parameters can not explain the deficits seen in the other experiments, so a hypothesis 
tests based on Xdn was performed to obtain an exclusion region, instead of a method based 
on a distance from the minimum \ 2 - Figures and |3] show shaded regions where the 
exclusion probabilities are more than 99% (x 2 > 15.09 for 5 dof), using the expected nadir 
dependencies for u e — > u^ T and u e — » v s , respectively. The typical parameters shown in 
figure [1] are excluded. Absolute flux information was not used in these calculations so the 
excluded regions found are independent of the 8 B neutrino flux. 

A combined rate analysis was performed using data from the 4 different solar neutrino 
experiments (CI, SAGE, Gallex and SK) with the assumption of the absolute fluxes given 
in ||. The analysis follows the method given in [JTJJ with updated theoretical uncertainties 
( 37 Cl cross section |16 |, 71 Ga cross section |17| and diffusion 0) and the latest experimental 



results [Q. In order to see the impact of the DN information obtained at SK, the rate 
analysis was performed with and without the DN information from SK. Results from this 
combined rate analysis are also shown in figures |2] and |[ Regions inside the solid and 
dotted lines show the allowed region at the 99% C.L. (% 2 < Xmin + 9.21) with and without 
DN information, respectively. 

In summary, SK observes no DN variation in the 8 P> solar neutrino flux in 504 days of 
data and sees no significant zenith angle variations. This result sets an exclusion region in 
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the oscillation parameter space where regeneration of v e would be expected. That exclusion 
region constrains the MSW solutions to larger mass difference (large angle solution) or to 
smaller mixing angle (small angle solution). 

We gratefully acknowledge the cooperation of the Kamioka Mining and Smelting Com- 
pany. This work was partly supported by the Japanese Ministry of Education, Science, 
Sports and Culture and the U.S. Department of Energy. 
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TABLES 



TABLE I. Day/Night fluxes obtained from 503.8 days of SK data. Flux values are normalized 
to 1 AU. Quoted systematic error for the "All" data is an absolute error but others are relative 
errors. 



Data set 


Nadir of the Sun 


Flux (10 6 / cm 2 /sec) 


Syst. error 


Live Time 


Day 
Night 
Nl 
N2 
N3 
N4 
N5 


-1 < cos(9 2 < 
< cos Z < 1 

< cos6 z < 0.2 
0.2 < cos6 z < 0.4 
0.4 < cos6 z < 0.6 
0.6 < cos6 z < 0.8 

0.8 < cos9 z < 1 


2.369™ (stat) 
2.481$™ (stat) 
2M0±%{$(stat) 
2.289ig-^(sta*) 
2.689±°ijt(stat) 
2.526+^ilKstat) 
2.3l8t° \H(stat) 


+U.6 % 
-0.5 /0 

-0.5 /o 
+1.3% 

+1.3 % 
-1.2 /0 
+1.3 % 
-1.2 /0 
+1.30/ 
-I.2/0 
+1.307' 
-1.2 /o 


242.2 days 
261.6 days 
43.6 days 
48.6 days 
64.3 days 

52.5 days 

52.6 days 


All 


-1 < cos6 z < 1 


2.436lffi(*tat)^ 


{syst) 


503.8 days 
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FIG. 1. Measured day/night solar neutrino fluxes as a function of the nadir of the Sun. Error 
bars represent statistical errors only. Night data is divided into 5 bins. Dotted histogram is the 
expected variation of a typical large angle solution and dashed histogram is that of a typical small 
angle solution. 
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FIG. 2. Flux independent exclusion region by SK day/night variation for v e — > z^ jT oscillations. 
Exclusion probabilities larger than 99% are shown in the shaded area. Regions inside of the dotted 
lines are allowed at the 99% C.L. from the combined rate analysis of Homestake, SAGE, Gallex and 
SK-flux in comparison with the BP98 SSM ||. Regions inside of the thick solid lines are allowed 
at the 99% C.L. from the combined rate analysis of the rates and the SK D/N variation. 
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